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Mechanical Assist Devices 
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• Motivation 


Over 5 million Americans and 20 million people worldwide suffer 
from Congestive Heart Failure (CHF) 

CHF pctients are still treated with drug therapy, 
however, at late stage heart transplantation is traditionally 
the only treatment hope 




Mechanical Assist Devices 
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• Motivation 




- Need for assist devices is very high 

Need : 25,000- 60,000 / YR 

Donor hearts available : 2,000-2,500 / YR 

(e.g. more than 4,000 patients were on the waiting list in 1999) 

- Need to find right match 

- Heart pump or VAD, for ventricular assist device, is being used 
as a temporary support to sick ventricle 

" BRI Di5E-TO-TRANSPLA NT" 


Mechanical Assist Devices 
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• Motivation 

- VAD vs Drug treatment, 

recent study suggests that 

• Survival rate for VAD patients vs for patients receiving drug treatment 
After 1 year 52 % VS 24.7% (it also depend; on the methods and drugs used) 
After 2 years 22.9 % vs 8.1 % 

• Some patients who stayed in ICU because of short of breath can walk a 
block after 1 year assisted by VAD 




Mechanical Assist Devices 
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• Motivation 

- VAD vs Drug treatment 

• However, complication rate for VAD is 2.35 times higher than that 
for drugs 

Comp icat ions include infections, bleeding, and mechanical malfunctions 
like motor failure, deformed tube and worn bearings 

Design improvements are needed to lower the risk, 

anc! possibly to use it as a permanent therapy (long-term device) 

''BRIDGE-TO-RECOVERY'' 


Mechanical Heart-Assist Devices 
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• Heart Volves 


• Ventricular Assist Device (VAD) 

Pulsatile Pump 

- Piston Driven : Low speed. Bulky 

- Pneumatically Driven : Need external support equipment 

Rotary Pimp 

- Axial F : low Pump : High speed. Small 
=* DeBctkey VAD is based on this concept 


• Total Artificial Heart 




Ventricular Assist Device 



• Requirements 

- Simplicity and Reliability 

- Small size for ease of implantation 

- Supply 5 liter/mi n of blood against 100 mmHg pressure 

- High pumping efficiency to minimize power requirements 

- Minimum Hemolysis and Thrombus Formation 





Computational Issues 
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• Geometry / grid definition 

Moving boundary 

• Solver 

Time accurate solver 

• Physical modeling 

Newtonian vs non-Newtonian 
Turbulence 

• Experimental & clinical data 


Solver: Viscous Incompressible Flow 
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• Formulation 

Can be viewed as a limiting case of compressible flow where the flow 
speed is nsignif icant compared to the speed of sound (Preconditioned 

compressible N-S eq.) 

=> Artificial compressibility approach 
- Artificial Compressibility Method (Chorin, 1967) 

INS3D family of codes 
Merkle et. d 
many more 

Or truly incompressible 

Pressure projection approach 

- MAC (Harlow and Welch, 1965) 

- Fractional Step Method (Chorin, 1968; Yanencko, 1971; Marchuck, 1975...-) 

- SIMPLE type Pressure Iteration (Caretto et d , 1972; Patanka A spalding, 1972...) 
=>Use der ived variables 

- Vorticity- Velocity (Fasel, 1976; Dennis et a!. ( 1979; Hafez et al., 1988 

- Stream function-vorticity 






Artificial Compressibility Method 
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• Formulation 



- Introduces hyperbolic behavior into pressure field. 

Speed of pressure wave depends on the artificial compressibility parameter , p. 

- The equations are to be marched in a time like fashion until the 
divergence of velocity converges to zero. 

=> Relaxes incompressibility requirement. 

Time variable during this process does not represent physical time step. 

For time-cccurate solutions 

- Iterate the equations in pseudo-time level for each time step until 
incompressibility condition is satisfied. 

=> Efficient sub-iteration is the key for success 


Artificial Compressibility Method (INS3D-UP) 
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• Time accurccy is achieved by subiteration 

- Discretize the time term in momentum equations using second-order 
three- point backward- difference formula 


3 q m " -4 q” + q* 
2 At 


= —(rhs)'*' 


- Introduce a pseudo-time level and artificial compressibility, 

- Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied. 
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• Code performance 

- Computing time : 50-120 ms/grid point/iteration (on C90 single cpu) 

- Memory usage: Line -relaxation 45 words/grid point 

GMRES-ILU(O) 220 words/grid point 
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A Pressure Projection Method(INS3D-FS) 
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• Approach in generalized coordinates 

- Finite volume discretization 

- Accurate treatment of geometric quantities 

- Dependent variables - pressure and volume fluxes 

- Implicr time integration 

- Fractional step procedure 

Solve auxiliary velocity field f irst, 

then erforce incompressibility condition by solving a Poisson 
equa-ion for pressure. 

• Code performance 

- Computing time ; 80 ms/grid point/iteration (on C90 single epu) 

- Memory usage: 70 words/grid point 


Pressure Projection Method 
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' Fractional-step 

- Solve for the auxiliary velocity field, using implicit predictor step- 

— (u -W) = -Vp" +h(u') 

At 

- The velocity field at time level (n+1) is obtained by using a 
correction step: 

— = -Vp~‘ + A(«r ) - Vp* + h(u‘ ) 

A t ' 

- The incompressibility condition is enforced by using a Poisson equation 
for pressure (p' = p'" -p') 

V 2 p' = — V-m' 

F At 
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^ Current Challenges 

r him Heseertfi Center _=======— ========= ** =s ^=* 1 ,=== !a! ===== — 

• Challenges where improvements are needed 

- Time-integration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 

• As the computing resources changed to parallel and distributed 
platforms, computer science aspects become important such as 

- Scalability (algorithmic & implementation) 

- Portability, transparent coding etc. 

• Computing resources 

- “Grid" computing will provide new computing resources for 
problem solving environment 

- High-fidelity flow analysis is likely to be performed using “super 
node' which is largely based on parallel architecture 
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• iNS3b-MPi 

(coarse grain) 

T. Faulkner <$ J. Dacles 


Group 1 Group 2 


• INS3D-MFT / Open MP 

MPI (coarse g-ain) + OpenMP (fine grain) 
Implemented jsing CAPO/CAPT tools 
H. Jin A C. Kir is 


• INS3D-MLP 

C. Kiris 



MPI 
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A OpenMP 
^■^^■^hreads 
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threads I 
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MLP Pfumw 2 
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Validation-SSME Turbopump Flow Analysis 
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• 5SME HPFTP 11' Impeller 

Shrouded impe ler: 6 full blades, 6 long partials, 12 short partials 6322 rpm, Re=1 ,81x10 s per inch 


HUB SURFACE COLORED BY STATIC PRESSURE 

Pressure 


COMPARISON WITH EXPERIMENTAL DATA 
IMPELLER EXIT PLANE AT 51% BLADE HEIG HT 

— * — Experiment I 
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Parallel Implementation of INS3D 



[NS3D-MLP (NAS MLP no pin-to-node) 
/ OpenMP 



TEST CASE : SSME Impeller- 
60 zones/ 19.2 Million points 
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Scripting for Acceleration of Grid Generation 
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INLET GUIDE VANES AND DIFFUSER 

Old IGV New IGV Old DIFF New DIFF 
No. of points (million) 7.1 1.1 8-0 1 - 6 

Time to build 1/2 day 10 sec. 1/2 day 8 sec. 

Script timings on new grids based on SGI Rl2k 3 00 MHz processor 
Time to build script = 1 day for IGV, 1 day for DIFF 




Parallel Implementation of INS3C) 
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INS3D-MLP / 40 Groups 


RLV 2 ml Gen Turbo pump 
114 Zones / 34.3 M grid points 


a 400 
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Numbor ol CPUs 


Per processor Mflop is 
between 60-70. 

Code optimization for 
cache based platforms 
is currently underway. 
Target Mflops is to 
reach 120 per 
processor. 

increasing number of 
OpenMP threads is also 
the main objective for 
this effort. 




SSME-rigl / Initial start 
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Time Step 13: Impeller rotated 8-degrees at 100% of design speed 



PRESSURE 


VELOCITY MAGNITUDE 


SSME-rigl / Initial start 


After 2 1/2 rotation: 


• Status 

- 34.3 Million Points 

- 400 physical time steps in one rotation. 

- One physical time-step requires less then 
12 minutes wall time with 128 CPU's on 
Origin platforms. One complete rotation 
requires 3.5-days wall-clock time with 
128 processors dedicated to the task. 

- I/O and memory management are critical 
for wall-dock time reduction 


Particle Traces and Pressure 


• Issues / Needs 

- In reolitv, more than 10% of the 
supercomputing facility to one task is not 
always practical. 

- Need lOOx bigger supernode or use lower- 
fidelity method 

- Communication to/from designers and 
experimental group is a part of critical 
technologies (in Orid computing) 








DeBakey VAD™ Pump Schematic 
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I INDUCER/IMPELLER 


FLOW STRAIGHTENER 
(FRONT BEARING SUPPORT) 


- MAGNETS EMBEDDED IN 
IMPELLER BLADES 


r DIFFUSER 

(REAR BEARING SUPPORT) 


VAD BEARING S ^ 

ZIRCONIA SHAFTS IN OLIVE 

RINGS, AGAINST SAPPHIRE ENDSTONES 


MOTOR STATOR 


Issues in Axial flow VAD 
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1 Problems Related to Fluid Dynamics 

- Small size requires high rotational speed 
Highly efficient pump design required 

-High shear regions in the pump may cause excessive blood 
cell damage 

Minimize high shear regions 

-Local regions of recirculation may cause blood clotting 
Good wall washing necessary 

iSmall size and delicate operating conditions make it difficult 
to quantify the flow characteristics experimentally 






NASA/DeBakey VAD (Baseline Design) 
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Flow Pattern Near Suction and Pressure Sides of Full Blade 

Tip Clearance : 0.009 in. 


Design I 


Traces Colored by Axial Velocity Magnitude 
-0.690 -0.365 -0.040 0.285 0.610 


Rotational Speed .12.600RPM 


Flow Rate : 5 lit/nun 


DeBakey VAD Development Timeline 
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• CFD /Assisted Design 

1993 - NAS A/ ARC is asked to develop CFD procedure to improve design and 
performance. D. Kwak and C. Kiris visit JSC to study the device 
The technology developed for rocket engine such as the Space Shuttle 
main engine was to be extended to blood flow Simulation 

1994 - Kir s and Kwak begin work on 
design analysis using NAS 
supercomputers BudirnDw 


, NEW DESIGN WAS PROPOSED TO 
INCLUDE AN INDUCER BETWEEN 
THE FLOW 5TRAIGHTNER AND 
THE IMPELLER 




Particle Traces Colored by Velocity Magnitude 
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DeBakey VAD - In situ 





DeBakey VAD Development Timeline 
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• Human Implantation in Europe 

1998 - On Movember 13, 1998, the first six DeBakey VAD s are implanted in 
European patients by Roland Hetzer and DeBakey at the German 
Heart Institute of Berlin. One of the patients, fifty six year old Josef 
Pristov, is able to return home and spend Christmas with his wife after 
a month's stay for recovery and monitoring at the clinic 




DeBakey VAD Development Timeline 
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• Human Implantation in USA 


1999 - US latent is granted for the device on September 9 , 1999 

January 2001 - 124 patients have received the device 

The longest successful trial period to date in human is 
approximately 1 year 
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NASA/DeBakey VAD 
Accomplishments to date (4/1/02) 


• 160+ pat ents implanted 

• Number of patients currently ongoing with device 



U5 trial 

Approved for trials involving 178 patients (14 male, 6 female have 
been performed) in a multi-center trial 

European trial 

Receivec "CE mark" (the EU equivalent to FDA approval) 

The VAD was implanted in 115 people with no device failure 

Results to date 

Favorable compared to existing VADs 

Small incidence of thrombus is being investigated 
=* Further computational support is necessary 




Summary and Discussion-1 
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Computational approach provides 

- a possibility of quantifying the flow characteristics: especially valuable 
for analyzing compact design with highly sensitive operating conditions 

- a tool for conceptual design and for design optimization 

CFD + rocket engine technology has been applied 

- to modify the design of NASA/DeBakey VAD which enabled human 
implantation 

Computing requirement is still large 

- Unsteady analysis of the entire system from natural heart to aorta 
involves several hundred revolutions of the impeller 

- During on heart beat, impeller has 125 revolutions 

- With 102 4 processors of Origin, one simulation (with several heart beat) 
from heort to aorta can be completed in months 

=> Received 2001 NASA Commercial Invention of the Year Award 
(6/25/2002) 


Summary and Discussion-2 
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Further study is needed 

- to assess long term impact of mechanical VAD on human body, which 
requires modeling flexible wall and non- Newtonian effect and 
better downstream boundary conditions 

There exist some gaps between 

- CFD (assuming IT is a part of CFD applications) and 
biomedical expertise 




